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Abstract 



Over times shorter than that required for relaxation of enthalpy, a liquid can exhibit 
striking heterogeneities. The picture of these heterogeneities is complex with transient 
patches of rigidity, irregular yet persistent, intersected by tendrils of mobile particles, 
flickering intermittently into new spatial patterns of motion and arrest. The study of 
these dynamic heterogeneities has, over the last 20 years, allowed us to characterize 
cooperative dynamics, to identify new strategies in controlling kinetics in glass-forming 
liquids and to begin to systematically explore the relationship between dynamics and 
structure that underpins the behaviour of amorphous materials. Computer simulations 
of the dynamics in atomic and molecular liquids have played a dominant role in all 
of this progress. While some may be uneasy about this reliance on modelling, it is 
unavoidable, given the amount of microscopic detail needed to characterize the dy- 
namic heterogeneities. The complexities revealed by these simulations have called for 
new conceptual tools. In this essay, I have tried to provide the reader with a clear and 
complete account of how these tools have been developed in terms of the literature on 
kinetic length scales in molecular dynamics simulations. Through the 'prism' of these 
length scales, this essay addresses the question what have we learnt about dynamic 
heterogeneities from computer simulations? 

0.1 Introduction 

What is it that distinguishes a glass from a crystal? Starting with the most casual 
inspection, the presence of oriented planes, grain boundaries or edges would indicate a 
crystal since the glass must be isotropic (and, hence, amorphous). Looking closer, the 
presence of sharp intensity peaks in scattered radiation at large angles indicates the 
presence of repeated parallel planes of density. The absence of such planes is typically 
all it takes for us to label the material 'disordered'. At the level of the constituent par- 
ticles (atoms or molecules), however, the most striking physical consequence of being 
in a glass rather than a crystal is the large number of different local environments (i.e. 
structural heterogeneity) in the former. Daama and Villars (Daama and Villars, 1997) 
have established that over 90% of the known 17000 inorganic crystal structures with 
intermetallic structure types consist of no more than 4 distinct coordination environ- 
ments. This is in clear contrast to a corresponding glass. Reverse Monte Carlo analysis 
of EXAFS measurements of the intermetallic glass N180P20 (Luo and Ma, 2008) have 
identified over 15 different coordination environments. Collective dynamics serves to 
amplify these structural variations, translating the often subtle differences in local con- 
figurations into dramatic variations in local relaxation kinetics. The result is that the 
approach to the glass transition is typically marked by striking spatial heterogeneity 
in dynamics. Many of the landmark insights of material science - dislocation-mediated 
plasticity, the Peierls barrier, Nabarro-Hcrring creep, grain boundary mobility - can 
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be comfortably re-expressed as manifestations of dynamic heterogeneities. In crystals, 
these objects are identifiable as defects and therefore have explicit structural signatures 
but, by shifting our focus to the spatial distribution of dynamics rather than order, we 
can expand our study of localised relaxation processes to materials for which we have 
no a priori notion of the structural origin of the localisation. From this perspective, 
dynamic heterogeneities provide a genuine opportunity to develop a universal descrip- 
tion of dynamic localization and collective relaxation in condensed matter - ordered 
and disordered, in equilibrium or out. 

Over the last 20 years, there has been a steady growth in the appreciation of the 
ubiquity of dynamic heterogeneities in disordered materials and of the value of their 
study. The subject of dynamic heterogeneities has been considered in a number of 
reviews of the glass transition (Poole, 1998; Kob, 1999; Glotzcr, 2000; Jackie, 2002; 
Andersen, 2005; Binder and Kob, 2005; Heuer, 2008). These spatial fluctuations, which 
take the form of transient kinetic domains, represent an extension of the 'traditional' 
phenomenology of disordered materials ('traditional' referring to thermodynamics and 
bulk averaged scattering and dynamic susceptibilities); an extension capable of pro- 
viding spatial information about the fluctuations associated with the collective dy- 
namics without requiring any insight as to the particle arrangements responsible. The 
existence of heterogeneities can also provide a link between different aspects of the 
phenomenology - stretched relaxations, fragility, etc. (Perera and Harrowell, 1996). 
While computational methods may dominate their study and theoretical goals pro- 
vide much of the motivation, it is worth emphasising that dynamic heterogeneities are 
not a theoretical construction but a physical fact and their description is, in the end, 
of value in its own right. 

To talk about dynamic heterogeneities in anything other than pictures we need 
some quantities to measure and kinetic lengths represent the most versatile of these, 
making them something of a lingua franca of glassy dynamics. What do we mean by 
a kinetic length? Prior to the study of dynamic heterogeneities, the kinetic length 
was attributed, somewhat vaguely, to a size of collective motion or a cooperative re- 
arrangement. In this essay we shall review what we now know about kinetic length 
scales associated with dynamic heterogeneities as a result of molecular dynamics (MD) 
simulations. The restriction to MD simulations, while omitting many important as- 
pects of the research into dynamic heterogeneities, ensures that we shall only consider 
those aspects of cooperative dynamics that we can explicitly connect to the positions 
and momenta of particles. We shall not cover the work on lattice models of glass 
forming liquids and we shall not cover the various theoretical treatments of the glass 
transition. Interested readers can find excellent reviews of these important topics in 
refs. (Ritort and Sollich, 2003) and (Cavagna, 2009), respectively, and in the chapters 
of this volume. 

We shall take the view that MD simulations confront us with all of the essential 
complexities of real supercooled liquids and glasses without, necessarily, providing a 
quantitatively accurate model of any one specific glass former. The great boon of 
these simulations, i.e. the provision of all particle positions and momenta at as many 
instances as desired, is also the source of the rather stringent obligation that they exert. 
After all, if you can access all information, then the failure to explain any aspect of 
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the glass transition that can be captured by a simulation can only be attributed to our 
personal failure in asking the right question. Dynamic heterogeneity owes much of its 
existence as a recognised phenomenon to MD simulation studies. While we now have 
mesoscopic (Weeks et ai, 2000; Abate and Durian, 2007) and macroscopic (Dauchot 
et ai, 2005) analogues of liquids in which the dynamic heterogeneities can be directly 
seen, it is fair to say that these discoveries owe their validation as models of microscopic 
heterogeneities to comparisons with computer modelling. 

The idea of a kinetic length scale, which extends back, at least, to the work of 
Adam and Gibbs (Adam and Gibbs, 1965) in 1965 (if not further back (Jenckel, 
1939)), is that the growing degree of dynamic correlation in a supercooled liquid can 
be measured in terms of a length that grows with cooling. To be accurate, Adam-Gibbs 
did not talk of a length scale but, rather, a number of particles required for a collective 
rearrangement. Converting this number into a length requires some sort of ancillary 
assumption concerning the connectivity of these particles. There is an undeniable 
wishfulncss involved in expecting to be able to replace something so poorly defined as 
the 'extent of cooperative motion' with something as concrete as a length. Computer 
simulations have provided the opportunity (obligation, really) to put this aspiration 
to the test and develop explicit expressions (theoretical or algorithmic) for this kinetic 
length scale. 

MD simulations require interparticle potentials. Simulations of supercooled liquids 
require interactions that are computationally simple in order to get as many time 
steps out of a computer as possible and yet geometrically complex enough to stave 
off crystallization long enough to allow characterization of the metastable liquid. In 
Table 0.1 we provide details of a number of models, all based on binary atomic alloys, 
that have proved popular and which form the basis of many of the simulation studies 
of the glass transition. 

In this review we shall consider five approaches to calculating a kinetic correlation 
length from MD simulations: direct measures of the spatial distribution of mobility, 
4-point correlation functions and the associated susceptibility X4i finite size effects, 
kinetic correlations at interfaces and the crossover lengths for kinetic properties. 

0.2 Kinetic Lengths From Displacement Distributions 

Historically, the first approach (Deng et ah, 1989) to analysing dynamic heterogeneities 
was to make a map of them. Maps are appealing. They retain a large amount of 
information. A single kinetic length scale, as we shall see, represents a major (and 
often uncontrolled) discarding of most of this information. The appeal of a map must 
be weighed against their principal short coming - they provide information (a lot of it, 
admittedly) about a single instance of heterogeneity only. Considerable care, therefore, 
is generally required in extracting the statistically significant features from the noise. 

The earliest example of displacement maps being used to specifically characterise 
dynamic heterogeneities can be found in a 1989 paper by Deng, Argon and Yip (Deng 
et ai, 1989). This paper contains a section with the prophetic title Inhomogeneities 
and the clustering of atomic motions. Inspection of maps of the particle displacements 
over a chosen time interval (see Fig. 0.1 for an example), reveal the generic features of 
the heterogeneities. Typically, we find connected domains of slow and of fast particles, 
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Label Inter-particle Potential Composition Reference 



BMLJ1 <Pij(r)=4e ij ([cri j /r} 12 -[ai j /r] 6 ) A 4 B a 

(Kob- Andersen) cjaa = 1-0, <?bb — 0.88, gab = 0.8 
6 A A = 1.0,e B B = 0.5,£AB = 1-5 



BMLJ2 <Pn(r)=4e ij ([ ( T ij /r] 12 -[a ij /r] 6 ) AB b 

(Wahnstrom) o~aa — 1-2,(Tbb = 1-0, (tab = 1-1 

e AA = tBB = tAB = 1-0 



SS Pij(r) =cij(cnj/r) 12 AB c 

fso/i sphere) gaa = l-2,a B B = 1-0, ctab = 1.1 

£AA = CBS = £AB = 1.0 



SD ^(r) =e ij (a 4 j/r) 12 (in 2D) AB d 

f.S0/£ dzsfcj (7AA = 1.4,CTBB = 1.0, (TAB = 1-2 

CAA = CBS = £AB = 1-0 



References: a) (Kob and Andersen, 1995), b) (Wahnstrom, 1991), c) (Bernu et at, 
1987), d) (Perera and Harrowell, 1998) 

Table 0.1 Details of a number of model glass forming liquids. Note that BMLJ1 and BMLJ2 
differ in that the former includes a strong preference for AB neighbours while the latter does 
not impose any specific chemical ordering. The soft sphere (SS) model is the same as the 
BMLJ2 system except for the absence of the attractions. 

with the former typically more compact than the latter. Collective strains within the 
'slow' domains coexist with low dimensional flows among the more mobile particles. 
Forced to identify lengths we might choose the average extent of slow domains or the 
fast domains or we might consider the length scale over which displacement direction 
is correlated. Alternatively, we might ask about the size of the 'core' regions where 
large displacements appear to be directed randomly or, instead, determine the average 
separation between such localised reorganizations. Our choice would, of course, be 
simplified if an argument existed that established that all of these lengths scaled in a 
similar fashion. Unfortunately, no such argument exists. 

In 1995, Hurley and Harrowell (Hurley and Harrowell, 1995) extracted a kinetic 
length scale for the one component soft disk liquid in terms of the decay in the variance 
of relaxation time maps as the linear dimension of the scaling volume increased using a 
box scaling method. (That these calculations were carried out on an equilibrium liquid 
- and only later on a supercooled mixture (Perera and Harrowell, 1998) - underscores 
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Fig. 0.1 The particle motions fi(t + 50r Q ) — fi(t) where r Q is the structural relaxation time 
for a polydisperse mixture of hard disks. [Reproduced with permission from ref. (Doliwa and 
Heuer, 2000). ] 

the point that dynamic heterogeneities are not restricted to supercooled liquids.) The 
relaxation time was defined as the first passage time for a particles displacement to 
exceed a threshold distance that was chosen to maximise the kinetic length. The kinetic 
length in the SD liquid was found to exhibit a super- Arrhenius increase with decreasing 
temperature (Perera and Harrowell, 1998; Perera, 1998; Perera and Harrowell, 1999). 

In 1997 Kob et al (Kob et ai, 1997) introduced an analysis of dynamic hetero- 
geneities based on determining the statistics of clusters of various kinetic subpopula- 
tions. Mobile particles were defined in such a way that they comprised of the ~ 5% 
of particles with the largest displacements over a time interval corresponding to the 
maximum in the non-Gaussian parameter (see Section 0.6). A characteristic size was 
obtained as S, the mass weighted average cluster size. For the BMLJ1 liquid, the tem- 
perature variation of S was fitted as S = 0.975/(T- 0.431) ' 687 (Donati et ai, 19996). 
Previously, the mode coupling theory (Gotze and Sjogren, 1992) had predicted the 
divergence of the relaxation time via a similar functional form, i.e. (T — T c )~^ . In the 
BMLJ1 liquid, T c ~ 0.43 (Kob et ai, 1997). The analogous average cluster size for the 
5% slowest particles exhibited little variation with T. The authors suggested that the 
apparent divergence of the size of the mobile clusters was significant, possibly linked 
to the divergence predicted by the mode coupling theory [32]. While an increase in 
the kinetic length with cooling is confirmed by all approaches, some care needs to be 
taken in interpreting an increase in cluster size when one is looking at clusters of some 
subpopulation of fixed size. A decrease in the total number of such clusters will force 
an increase in the size of those remaining, simply as an artefact of how the clusters 
are defined. 
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The cluster analysis has been applied to a number of systems: polymers (Bcmie- 
mann et at, 1999; Gebremichael et ai, 2001), water (Mazza et at, 2006), SiC>2 (Vogel 
et ai, 2004; Vogcl and Glotzer, 2004), the Dzugutov potential (Gebremichael et ai, 
2004) and polar diatomic molecules (Palomar and Sesc, 2008). Giovambattista et 
al (Giovambattista et ai, 2005) have calculated the fractal dimension of the mobile 
clusters in supercooled water. They found that large clusters exhibited a fractal di- 
mension of ~ 2, a value similar to that predicted (Lamarcq et ai, 2002) for low 
energy excitations in a spin glass. Vollmayr-Lee et al (Vollmayr-Lee et ai, 2002) have 
extended the mobile cluster analysis to a BMLJ1 mixture below its glass transition 
temperature. Below T g , mobility is determined using the mean amplitude of a parti- 
cle's fluctuation about its mean position. The authors report that clusters of mobile 
particles defined by this means were more compact and shorter lived than the anal- 
ogous clusters above T 9 . As to the question of the connection between the kinetic 
lengths of mobile particles on either side of the glass transition, strong correlations 
have been demonstrated (Widmer- Cooper and Harrowell, 2006) above T g in the SD 
mixture between the spatial distribution of mobile particles and those particles which, 
over short times, exhibited large amplitude fluctuations (analysed in terms of local 
Debyc- Waller factors (Widmcr-Coopcr and Harrowell, 2006)) - a criterion similar in 
spirit to that used by Vollmayr-Lee et al (Vollmayr-Lee et ai, 2002). 

Having introduced the cluster analysis, Donati et al (Donati et al, 1998) aug- 
mented it in order to examine the correlation between displacement directions and the 
positions of mobile neighbours. Starting with the same fraction of mobile particles, 
they introduced an additional requirement for belonging to a cluster - for two neigh- 
bouring mobile particles to belong to the same cluster the new position of one particle 
has to lie within a selected distance of the neighbours old position. This new overlap 
condition strengthens the interpretation of the clusters as representing a cooperative 
motion rather than simply some general aggregation. Striking pictures of string-like ar- 
rangements of displacements (i.e. collective movements that are correlated only along 
the direction of flow) have been presented for the BMLJ1 (Donati et ai, 1998) and 
SS (Kim and Yamamoto, 2000) mixtures. The observation of low dimensional particle 
flows raises some interesting questions about the mechanisms of cooperative dynamics 
and readers are encouraged to read refs. (Vogel et ai, 2004; Gebremichael et ai, 2004) 
closely for details of the complex motions by which displacement strings are formed. 

Since first presented in 1998 (Donati et ai, 1998), images of strings of particle dis- 
placements have proved a popular leif-motif for dynamic heterogeneities in general. 
To what degree is this representation accurate? Are dynamic heterogeneities gener- 
ally string-like or are extended strings of displacements simply the more eye catching 
members of a broad distribution of cluster shapes? There is not a lot of information 
that addresses this specific question. Focusing on the most mobile ~ 5% of particles 
in a BMLJ1 mixture, Donati et al (Donati et ai, 1998) demonstrated that the dis- 
placement vectors exhibited a local forward alignment and that the clusters defined 
using the overlap condition were described by an exponential distribution of the size 
with the average number increasing from ~ 1.4 (T = 0.55) to ~ 2.2 (T = 0.45). The 
mass weighted average size is larger, reaching ~ 15 at T = 0.45 (Donati et ai, 1999&). 
While the authors of ref. (Donati et ai, 1999 b; Donati et ai, 1998) refer to these clus- 



Kinetic Lengths From Displacement Distributions xi 



ters as 'strings', the overlap condition for cluster membership establishes only a flow 
and does not specifically establish the dimensionality of that flow. While the overlap 
condition will certainly include string-like correlations, it will also include some col- 
lective flows and strains in general. The existence of such collective flows in liquids at 
rest has been known for some time. Alder and Wainwright (Alder and Wainwright, 
1967) demonstrated in 1967 the local forward alignment between particle velocities 
in an equilibrium liquid, part of a pattern that strongly resembled the solution of 
the Navier-Stokes equations about a moving particle. Doliwa and Heuer (Doliwa and 
Heuer, 2000) have reproduced this pattern in a study of the displacements in a binary 
hard disk mixture. To establish string-like correlations one would like to establish the 
average coordination number of particles in the mobile cluster is < 2. This quantity 
has not been determined but Donati et al (Donati et al., 19996) have reported a frac- 
tal dimension of 1.75 for mobile clusters in the BMLJ1 model (obtained over a single 
order of magnitude data set). There are indisputably string-like objects such as the 
self avoiding random walk in 3D with a similar fractal dimension (1.66 from Flory's 
argument (Flory, 1969)). There are, however, other random objects with similar frac- 
tal dimensions - on a 3D lattice we have the random cluster below percolation with a 
fractal dimension of 2 and the backbone of the percolating cluster with a dimension ~ 
1.7 (Stauffer and Aharony, 1994) - for which the description 'string-like' really doesn't 
apply. The cumulative data clearly points to mobile clusters that are not compact 3D 
objects but it leaves open the question of what they are. There is certainly a need for 
more systematic studies of the shape (and its dependence on temperature and the time 
interval used to define mobility) as well as the length scale of dynamic heterogeneities. 

In simulations of the SD and SS mixtures, Yamamoto and Onuki (Yamamoto and 
Onuki, 1997; Yamamoto and Onuki, 1998a; Yamamoto and Onuki, 19986; Yamamoto 
and Onuki, 1999) shifted the focus from particle displacements to changes in the 
topology of a configuration in the form of 'broken bonds'. Their analysis involved 
generating spatial maps of the positions where initially nearest neighbour pairs first 
moved far enough apart for the two particle 'bond to be considered broken. These 
authors showed that the structure factor S&(q) of the broken bonds in both 2D and 
3D, accumulated over a fixed time interval, obeyed the small q expansion of S^ 1 (q), 
i.e. 



i + £ <r 



where £ represents the associated length. (The functional form in Eq. 0.1 corresponds 
to an exponential decaying correlation between broken bonds in real space.) Observing 
that Sb(q) appeared independent of temperature (and, hence, £) at large q, the authors 
concluded that S&(0) ~ £ 2 , the Ornstein-Zernicke result, and noted the analogy with 
critical fluctuations. 

Analogous correlation functions have also been developed to describe the spatial 
correlations between displacements. In 1998, Poole et al (Poole et al., 1998) introduced 
the displacement-displacement correlation function G u (f,At), defined as follows. For 
the displacement amplitude field u(r,t,At) given by 
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N 

u(r,t, At) = J2 ^(t + At ) - ?i(t)W- n{t)) (0.2) 

i=l 

we can write 

G u (f,At)= J d?< [u(?+r,t, At)- < u >][u(f, t,At)- < u >] > (0.3) 

By analogy with the relationship between the density-density correlation function, 
the variance of the number of particle N and the compressibility k, Donati ct al (Donati 
et al., 1999a) write 

J dfG u (r, At) =< [U- < U >] 2 >=< u >< U > kTn u (0.4) 

where U = J dru{r, t, At) is the total displacement and k u is a time dependent dy- 
namic susceptibility. For the BMLJ1 model, the maximum value of k u with respect to 
the displacement time interval At varies with T as (T — 0.435)" 84 . From the decay 
of the density correlations with separations or, equivalently, working with the analo- 
gous structure factor and Eq. 0.1 as in ref. (Yamamoto and Onuki, 1998a), a kinetic 
length can be extracted. This kinetic length corresponds to the correlation length of 
fluctuations in displacement amplitude field. A number of models have been analysed 
using this approach: the BMLJ1 mixture (Poole et al. 1998; Donati et al., 1999a), 
molecular (Qian et al., 1999) and polymeric (Bcnnemann et al., 1999; Gebremichael 
et al., 2001) liquids and polydisperse hard spheres (Doliwa and Heuer, 2000). 

The connection between length scales and relaxation times represents the main 
motivation for looking at the length scales in the first place. Empirical power law 
relations between a relaxation time and the length of the form 

t = Ai z (0.5) 

were reported in all of the above cases, often with large exponents. In the two dimen- 
sional SD mixture z ~ 4 from both the box scaling (Perera and Harrowell, 1998) and 
broken bond structure factor (Yamamoto and Onuki, 1998a). An even larger exponent 
is found in the BMLJ1 mixture (Donati et al., 19996), where t q ~ S 4 ' 5 . Since S is 
the average (mobile) cluster size, it will be related to a length £ via S ~ where 
1 < v < 3, implying, when substituted back into the power law relation with r a , a 
value of z > 4.5. The BMLJ2 and SS mixtures, in contrast, have modest exponents: 
z ~ 2 (bond breaking) (Yamamoto and Onuki, 1998a) for the SS mixture and z ~ 2.34 
(4-point correlation) (Lacevic et al., 2003) in the BMLJ2 model. We remind the reader 
that the BMLJ2 and SS models have the same short range repulsions and associated 
length scales. The smaller the value of z, the larger the kinetic length required to 
achieve a given relaxation time. Does this mean that the structures responsible for £ 
in these two liquids are mechanically less robust than in the SD and BMLJ1 liquids, 
since they must be larger to achieve an equivalent stability? 

The difference between the value of z for the SS and SD liquids is also noteworthy 
since both only make use of short range repulsions. They differ, of course, in their 
spatial dimension. Yamamoto and Onuki (Yamamoto and Onuki, 1998a) comment 
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Method 


Selected 

length/time 

scales 


Manner of Assigning Values 


box scaling 


threshold length 


maximise kinetic length 


mobile cluster 6 


observation time 


maximise non-Gaussian parame- 
ter 

maximum length r* for which 

G(r*,t) =G gaus (r*,t) 




threshold length 


bond breaking c 


observation time 

maximum bond 
length 


0.5Tb where r& is the average 
bond breaking time 
a length lying between the first 
two peaks of gij{r) 


displacement-displacement 
correlation function" 1 


observation time 


maximise susceptibility K u (t) 



a (Hurley and Harrowell, 1995), b (Kob et al, 1997), c (Yamamoto and Onuki, 
1998a), d (Poole et al, 1998) 



Table 0.2 A summary of the various methods of calculating kinetic lengths and the associ- 
ated length and time scales required in each case. 

that the difference in z represents the only significant difference they observed in the 
features of dynamic heterogeneities between 2D and 3D. Where direct comparisons 
have been carried out (Doliwa and Heuer, 2000; Yamamoto and Onuki, 1998a), there 
is no evidence that the long-wavelength anomalies that are a well documented feature 
of crystals and liquids in 2D have any significant impact on the glass transition in 
2D. In this context, it is worth noting that the phenomenology of supercooled liquids 
in 3D and 4D is also very similar, with the most significant difference being that 
the breakdown of the scaling between the diffusion constant and the relaxation time 
(the Stokes-Einstcin relation) is somewhat weaker in the higher dimension (Eaves and 
Reichman, 2009). 

The kinetic length scales described in this Section - direct measures of the spa- 
tial distribution of particle displacements or bonds breaking - have proven valuable 
descriptive tools in establishing the reality of dynamic heterogeneities and their de- 
pendence on temperature. In Table 0.2 we provide a summary of these approaches, 
identifying the time scales and/or lengths that must be chosen to resolve the transient 
heterogeneities in each case and the methods used to assign values to these quantities. 
In an important development, explicit in box scaling and the displacement suscepti- 
bility approaches and implicit in the mobile cluster method, it becomes legitimate to 
select a quantity, not based on some particular physical justification, but rather on 
the pragmatic goal to maximise the resolution of the heterogeneity. 

Is there a best measure of a kinetic length? There does appear to be a consensus 
concerning the choice of how to obtain a kinetic length. The kinetic length of choice 
is that associated with %4, based on 4-point correlations and closely related to the 
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displacement-displacement correlations introduced by Poole et al (Poole et al, 1998). 
As we shall see in the following Section, xt has a number of appealing features, not 
least being its close connection with the formalism of spin glasses. Comforting as such 
consensus can be, it is worth pointing out that Xi provides the same information as 
that contained in the other measures reviewed in this Section. Its popularity depends, 
not on the superiority of its description of dynamic heterogeneities, but on its accessi- 
bility, particularly from experiments, and on its potential role in developing theoretical 
treatments. Whether X4 provides a sufficiently complete account of cooperative dy- 
namics is a question we shall return to in Section 0.7. 

0.3 Kinetic Lengths From 4-Point Correlations Functions 

A minimal description of dynamic heterogeneities requires that we measure the statis- 
tical correlations between the movement of pairs of particles. This description requires 
a correlation involving four positions: fi(t), fi(t + r), ^(i) and ri(i + r). The first 
calculation of such 4-point correlation functions from MD simulations of a Lennard- 
Joncs mixture (similar to BMLJ2 except the size ratio (y\i/<J22 = 5/8) was reported by 
Dasgupta et al (Dasgupta et al, 1991) in 1991. The motivation of these calculations 
was to test for the presence of a growing length scale associated with fluctuations 
of the Edwards- Anderson order parameter, lim t _>. 00 < 8n(f,t )5n(r,t + t) >, where 
n(f, t) is the density field. Such a growth in the length scale of an order parameter 
fluctuation would have provided evidence of a thermodynamic glass transition. Fixing 
the spatial separation (i.e. \r\ — fal) at 2cth, the authors found no evidence of corre- 
lated fluctuations at any time. It is likely that this negative result was a consequence 
of the choice of time interval. As we shall see, fluctuations of the density autocorrela- 
tion function exhibit a maximum at around r a , while the calculations of Dasgupta et 
al. (Dasgupta et al, 1991) were, by design, carried out in the plateau interval of the 
relaxation function, well short of this time scale. 

In 1999, Franz, Donati, Parisi and Glotzer (Franz et al., 1999; Donati et al., 2002) 
described how the displacement-displacement correlations introduced to quantify the 
spatial distribution of particle movement (see previous Section) could be reformulated 
in terms of 4-point correlations 1 This reformulation involved a conceptual convergence 
of a method introduced to describe the spatial correlations of displacements with the 
formal tools developed to treat fluctuations in spin glasses, such as motivated Dasgupta 
et al (Dasgupta et at, 1991). Beginning in 2000, Glotzer and coworkers (Lacevic et al, 
2003; Glotzer et al, 2000; Lacevic et al, 2002) presented a detailed account of the 
4-point correlation function formalism as applied to supercooled liquids. Here we shall 
summarise their expressions (Lacevic et al, 2003) for a kinetic length scale using a 
4-point correlation. 

The structural relaxation is described by the quantity 

N N 

et^EEH^c)-^')!) (°-6) 
»=i j=i 

1 To get the chronology straight during this busy period, note that the paper (Donati et al., 2002) 
first appeared as a preprint cond- mat/9905433 in 1999. 
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a measure of the degree to which a configuration at time t still overlaps the initial 
arrangement. Overlap of a particle with itself or another particle at an earlier time is 
established through the introduction of a window function w(r) (where w(r) = 1 if 
|r| < a and zero, otherwise). The self part Q s (t) of the relation function can be dehncd 

as 

TV 

Q s (t)=$>(|rH0)-f^)|) (0.7) 

i=l 

Lacevic et al (Lacevic et at, 2003) showed that this self part is the dominant con- 
tributor to the relaxation function, the dynamic susceptibility and the kinetic length. 
Physically, this means that the essential relaxation event corresponds to the departure 
of a particle from its own initial 'site'. Eq. 0.7 provides an interesting link with the 
previous approaches (Hurley and Harrowell, 1995; Kob et al., 1997) in which the de- 
scription of dynamic heterogeneities involved considering each particle moving beyond 
some threshold distance from its initial position. Similarly, Stein and Andersen (Stein 
and Andersen, 2008; Stein, 2007) examine the 4 point correlations of a mobility defined 
as 

»(r,t) = J2t(r-n(t)hi(t) (0-8) 

i=l 

where Hi(t) = 1 if |fj(t + 1*) — r(t)\ > d and zero otherwise. 

Toninelli et al (Toninelli et al., 2005) and Flenner and Szamel (Flenner and Sza- 
mel, 2007), along with others, have used the (microscopic) self intermediate scattering 
function F s (q, t) instead of the self overlap function Q s (t). There is no fundamental dif- 
ference between the two relaxation functions. The choice of the scattering wavevector q 
in F s (q, t) selects the reference length scale analogous to the choice of the value of a in 
the overlap function. The authors in ref. (Toninelli et al, 2005) considered how various 
models of collective behaviour (elastic modes, domain wall fluctuations, etc.) where 
represented at the level of the 4 point correlations. Flenner and Szamel (Flenner and 
Szamel, 2007) showed that the correlations of the fluctuations in the microscopic self 
intermediate function exhibited an anisotropy associated with the direction of particle 
motion relative to the vector between particle pairs. This latter resulted expressed, 
in terms of the 4 point correlation function, the anisotropy that had previously been 
established by Doliwa and Heuer (Doliwa and Heuer, 2000) who demonstrated that 
the kinetic length scale (obtained using Eq. 0.3) along the direction of the particle dis- 
placement was roughly twice as long (in a binary hard sphere mixture) as that along 
a direction perpendicular to the particle displacement. 

A susceptibility Xi(t)y analogous to n u (t) defined previously for the displacement 
correlations, is defined in terms of the fluctuations of the relaxation function Q(t), i.e. 

Xi{t) = ^{<Q 2 {t)>-<Q{t)> 2 ) (0.9) 

The expression for the susceptibility Xi(t) (Eq- 0.9) in terms of the variance of the 
structural relaxation function Q(t) is remarkable in that it provides information about 
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the extent of dynamic heterogeneities (see below) without ever requiring that the 
dynamics be spatially resolved (as the methods in Section 0.2 all did). Qualitative 
considerations of the fluctuations of Q(t) also provide a useful way of differentiating 
the order parameter fluctuations, envisioned by Dasgupta et al (Dasgupta et al, 1991), 
from the fluctuations in dynamics that are the subject of this review. In Fig. 0.2 we 
present sketches of two distinct types of fluctuations in Q. The top panel represents 
fluctuations of the height of the plateau. The Edwards- Anderson order parameter 
for spin glasses is this plateau height in an arrested system. An alternative type of 
fluctuation involves variations of the relaxation time (Fig. 0.2, middle panel). These 
two types of fluctuations are easily distinguished in their respective susceptibilities 
(Fig. 0.2, bottom panel). Where the plateau fluctuations result in a low amplitude 
X4(t) 7 extended over the whole plateau time region, the dynamic fluctuations typically 
produce a more sharply peaked X4(i), with the maximum occurring at roughly the 
structural relaxation time. Compare these two (idealized) possibilities with the Xi(t) 
calculated from simulations of the BMLJ2 mixture (Laccvic et al, 2003) in Fig. 0.3 
and wc find that fluctuations in dynamics, sketched as option b) in Fig. 0.2, pro- 
vide a reasonable correspondence with the simulated liquid. The implication is that 
it is the dynamic fluctuations, as opposed to those of the plateau/order parameter, 
that dominate the observed susceptibility in the supercooled liquid. Kirkpatrick and 
Thirumalai (Kirkpatrick and Thirumalai, 1988) were the first to point out that the 
fluctuations contributing to Xi(t) could arise from these two distinct sources: order 
parameter fluctuations and fluctuations in the dynamics themselves. 

The susceptibility Xi{t) can be directly obtained (Lacevic et al, 2003) from the 
4-point correlation function via 

X 4(t) = P j dr 9i (r,t) (0.10) 

where 



94(r,t) = ±(j2s(r-f k (o) + mM\m-m\M\Mo)-m\)) - 

P ijkl 

(0.11) 

or 

gi[ r,t)^gf{r,t)-(^j (0.12) 

where we have assumed that the correlations arc isotropic and so only depend on the 
magnitude of the separation. The significance of gf (r, t) is that it corresponds to pair 
correlation function of those particles at t = that end up overlapping with a particle 
at the later time t. It is this quantity, or rather its Fourier transform, 



S?{q,t)= / dfg?(r,t)exp(-iq-r?) 



(0.13) 
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Fig. 0.2 Sketches of the fluctuations in Q(t) due to a) fluctuations in the plateau height 
and b) fluctuations in the magnitude of the structural relaxation time. The time dependent 
dynamic susceptibilities Xiit), calculated using Eq. 0.7 for both types of fluctuations, are 
plotted in the bottom panel. 

that Lacevic et al (Lacevic et al., 2003) use to obtain the kinetic length £4 whose 
dependence on t is plotted in Fig. 0.4. As emphasised in ref. (Lacevic et al, 2003), 
considerable care needs to be taken in ensuring that the 4-point correlator used to 
extract the kinetic length does not include a weak 0(1/N) tail associated with bulk 
fluctuations. As is generally the case, fluctuations, such as measured by the 4-point 
correlations, are very dependent on the choice of ensemble. Dalle-Fcrricr et al (Dallc- 
Ferrier et al, 2007) provide a useful discussion of this point. 

Toninelli et al (Toninclli et al, 2005) have presented a thorough discussion of the 
behaviour of X&{fy over various time domains. The kinetic length increases in time as 
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Fig. 0.3 The time and temperature dependence of X4(i) for the BMLJ2 mixture. The time 
corresponding to the maximum in x*($) ls found to be similar in value to the structural 
relaxation time r a and to exhibit a similar T dependence. [Reprinted with permission from 
ref. (Lacevic et at, 2003) Copyright 2003, American Institute of Physics.] 



Fig. 0.4 The time and temperature dependence of the kinetic length ^(i) in the BMLJ2 
mixture. [Reprinted with permission from ref. (Lacevic et al, 2003). Copyright 2003, Ameri- 
can Institute of Physics.] 

the heterogeneous character of the unrclaxed domains is exposed by the dynamics until 
it reaches a maximum, beyond which is decays to zero as the persistent domains finally 
succumb to relaxation and homogeneity is re-established. The temperature dependence 
of the maximum length £ 4 (t*) in the BMLJ2 liquid is fitted by £ 4 (t*) = 0.82(T/T C - 
1)~ - 82 with an apparent divergence at the mode coupling temperature T c 0.57 (Lacevic 
et al., 2003). Such apparent divergences generally do not actually eventuate and there 
is a considerable literature discussing the various interpretations of both the power 
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law temperature dependence and the 'avoidance' of the singularity (Cavagna, 2009). 
As already mentioned, a scaling law, t q ~ £ 4 (i*) z , was found to hold with z = 2.34. 

The length scale a imposed by the window function w(r) plays a significant role 
in selecting the physical character of the fluctuations being measured. As pointed out 
in ref. (Lacevic et al., 2003), select a too small and the results are dominated by vi- 
brational motions that obscure any heterogeneities, select a too large and the notion 
of overlap quickly loses any physical significance as each particle can overlap with 
multiple particles. Dauchot ct al (Dauchot et al., 2005) have chosen a value for a the 
same way that t* is chosen, i.e. by finding the value that maximises Xi- Chandler ct 
al (Chandler et al., 2006) pointed out that the use of a small a (roughly, a/cr < 0.3) 
results in the spatial clustering of immobile particles dominating the observed fluc- 
tuations, whereas the use of a large a (i.e. 0.5 < a/a) results in a Xi(t) reflecting 
correlations in mobile particles. This length dependent crossover is related (Chandler 
et al., 2006) to the non-Fickian-to-Fickian crossover discussed in Section 0.6. Char- 
bonneau and Rcichman (Charbonncau and Reichman, 2007) have described how the 
dependence of Xi on the reference length scale differs when comparing liquids whose 
arrest is dominated by short-range repulsions and those dominated by short-ranged 
attractions. 

The great attraction of the 4-point correlation function approach to the kinetic 
length is not the 4-point correlation functions, g^r, t) or S±{q, t), which are at least as 
difficult to use to calculate a kinetic length as any of the other methods described in the 
previous Section. The real appeal is in the quantity X^(t) an( i related susceptibilities. 
As the space integral of the 4-point correlation function (Eq. 0.8), X4(**)i the maximum 
susceptibility, represents a 'volume of correlation' or, alternatively, a number N corr of 
correlated particles (Dalle-Ferrier et al., 2007). As the variance of the relaxation func- 
tion Q(t) (Eq. 0.9), Xi{t) is n0 more difficult to calculate than the relaxation function 
itself. This latter virtue marks the superiority of Xi(t) over the analogous susceptibil- 
ity, n u (t), from the displacement correlations (defined in Eq. 0.4). Comparing Figs. 0.3 
and 0.4, it is evident that x&(t) exhibits a qualitatively similar time dependence to 
that of the kinetic length £i(i) with a peak at some intermediate time, corresponding 
to a maximum in the differentiation of overlapping and non-overlapping domains. 

This theoretical accessibility has been extended towards experimental accessibility 
in a series of papers (Dalle-Ferrier et al., 2007; Berthier et at, 2005a; Berthier et al., 
2007a; Berthier et al., 20076; Biroli et al., 2006) starting with Berthier et al (Berthier 
et al., 2005a) in 2005. In ref. (Berthier et al., 2005a), the authors showed how a lower 
bound on the value of Xi(t*) could be obtained from 3-point correlations defined as the 
response of the structural relaxation function (i.e. Q(t) or its analogue) to a change in a 
control parameter, such as temperature, pressure or density. This connection between 
a 3-point correlations and N corr has been explored in some detail (Dalle-Ferrier et al., 
2007). Starting with the maximum of the 3-point susceptibility xt(£*) associated with 
the dynamic heterogeneities that are correlated to local enthalpy fluctuations, the 
number of correlated molecules is given (in the NPT ensemble) by 




(0.14) 
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Eq. 0.14 has been used to determine the size and temperature dependence of N corr 
for a range of molecular glass formers (Dalle-Ferrier et ah, 2007). The discussion in 
ref. (Dalle-Ferrier et ah, 2007) of these results raises a number of important questions 
regarding kinetic lengths in general. The correlation volumes at T g were found to be 
modest in size. This result offers hope that simulations can contribute to the descrip- 
tion of cooperative dynamics closer to T g than previously thought. It also cautions 
against arguments based on large separation of length scales in the supercooled liquid. 
Dalle-Ferrier et al (Dalle-Ferrier et ah, 2007) cast doubt on any simple connection 
between N corr and the cooperatively rearranging regions as imagined by Adam and 
Gibbs (Adam and Gibbs, 1965). There is a quite basic difficulty in trying to translate 
observed correlations in mobility into the mechanism responsible for those correlations. 
We shall return to this point at the end of the review. 

If both Xt(t*) an d £,4(t*) can provide information about the extent of dynamic 
heterogeneities, what, exactly, is the connection between them? If one assumes that 
the spatial distribution of the mobile particles is scaled by a single length in 
this case), it follows (Biroli et ai, 2006; Biroli and Bouchaud, 2004; Whitelam et ah, 
2004) that X4(t*) an d £,4(t*) are related by a power law, 

X4 (o = A^un 2 - 11 (o.i5) 

Stein and Andersen (Stein and Andersen, 2008; Stein, 2007) have confirmed the power 
law relation for the BMLJ1 mixture and found t] = —2.2. More recently, Karmakar 
et al (Karmakar et ah, 2010a) have carried out similar calculations but using a larger 
number of particles and reported a significantly smaller exponent, rj = —0.4. Karmakar 
et al (Karmakar et ah, 2010&) argue that the extraction of the kinetic correlation length 
from the 4-point structure factor Si(q,t) by fitting an Ornstein-Zcrnickc expression, 
i.e. 6*4(0, t) — Xo/(l + Q 2 ^ 2 ), is inaccurate unless large systems (i.e. N ~ 10 5 ) are used. 

Robust in definition and simple to apply, the susceptibility Xi{t*) has proved a 
popular measure of the extent of cooperative motion. As obtained from the fluctuation 
of the 2-point correlation, X4(t*)'s have been reported for colloids (Duri and Cipellctti, 
2006; Cipelletti et ah, 2003; Brambilla et ah, 2009; Ballesta et ah, 2008), granular 
material (Dauchot et ah, 2005; Keys et ah, 2007) and foams (Mayer et ah, 2004). 
In simulations, X^(t*) is being used to study cooperative behaviour in a wide range 
of systems - including those in non-Euclidean spaces and out of equilibrium. Some 
examples of the latter applications: Sausset and Tarjus (Sausset and Tarjus, 2010) 
have calculated Xi(t*) f° r a liquid of Lennard- Jones disks on a hyperbolic surface, 
Furukawa et al (Furukawa et ah, 2009) have analysed the 4-point correlation functions 
in a liquid under steady shear (finding mobile regions to be elongated) , and Parsaeian 
and Castillo (Parsaeian and Castillo, 2008) have studied the effects of aging on Xi(t*)- 
Abraham and Bagchi (Abraham and Bagchi, 2008) have demonstrated that the low 
temperature magnitude of Xi(t*) i n a polydispcrse Lennard- Jones mixture decreases 
with increasing width of the particle size distribution. 

0.4 Kinetic Lengths From Finite Size Analysis 

If heterogeneities are characterised by a length, then it follows that, as the size of a 
finite system approaches this inherent length, properties of the dynamics should show 
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a dependence on the system size. For this finite size effect to be linked specifically 
to a kinetic length scale it is important to also establish that the static correlations 
in the liquid are independent of the system size over the size range studied. In 1992, 
Dasgupta and Ramaswamy (Dasgupta and Ramaswamy, 1992) reported the absence 
of any size dependence in the relaxation time of a density autocorrelation function 
for a binary Lennard- Jones mixture with a radius ratio of 5/8. This study, however, 
failed to collect data for temperatures between 0.8T S and 0.99T g , despite what appears 
to be an onset of a finite size dependence in this temperature interval. The authors, 
expecting to see the time scale decrease with system size, interpreted the (slight) 
signs of the opposite trend as supporting their conclusion that there was no system 
size dependence. In 2000, Kim and Yamamoto (Kim and Yamamoto, 2000), studying 
the SS mixture for sizes N = 108, 1000 and 10000, found that as the temperature 
was lowered the relaxation time r Q exhibited a dependence on the system size, the 
relaxation time increasing as the system size decreased (see Fig. 0.5). No dependence 
of the pair distribution function on system size was found for all temperatures studied. 

The influence of system size on the relaxation times in the BMLJ1 mixture has been 
the subject of a number of studies. In 2003 Doliwa and Hcuer (Doliwa and Heuer, 2003) 
calculated the diffusion constants for four system sizes - N = 65, 130, 260 and 1000 - 
down to T c (~ 0.43). At the lowest temperature, they found D e5 /Di 3n ~ 1.2 but with 
a large error (±0.2). They concluded that the size dependence was small, significantly 
less (and in the opposite direction) to that observed in the SS mixture (Kim and 
Yamamoto, 2000). Stein (Stein, 2007), working with larger systems (N = 1000 and 
8000), found a slight decrease in the diffusion constant for the smaller system but, 
again, the difference was within the noise. In 2009, Karmakara et al (Karmakar et al., 
2009) determined the relaxation time in a large number of systems across the size 
range 50 < N < 1600 range of system sizes down to T = 0.45. They find a systematic 
increase in relaxation time with decreasing system sizes for N (see Fig. 0.9). The 
onset value of N below which this finite size effect is observed increases from ~ 100 
at T = 0.8 to ~ 200 at T = 0.45. These workers also studied the size dependence 
of Xi(~b*) and, using a method first adapted from the study of critical phenomena to 
the glass problem by Berthier (Berthier, 2003), determined a length scale from this 
data that increased on cooling from 2.1 (T=0.70) to 6.2 (T=0.45). An increase in the 
relaxation time with decreasing system size has also been reported in simulations of 
silica (Zhang et al., 2004; Teboul, 2006) but, unfortunately, there is no confirmation 
that the static properties remained unchanged, a nontrivial condition given the long 
range character of the potential in this model. 

With the exception of rcf. (Doliwa and Heuer, 2003), the data from simulations 
presents a picture of a modest but systematic increase in the volume associated with 
cooperative motion on cooling. Significantly, relaxation is slower in small systems. 
We shall consider the implications and possible origin of this behaviour in Section 
0.7. Confinement, whether imposed by walls or arising from inherent fluctuations in 
the supercooled liquid, appears to represent one generic mechanism for slowing down 
relaxation. 
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Fig. 0.5 The temperature dependence of r a for iV = 108 (open diamonds), 10 3 (closed 
diamonds), and 10 4 (open squares) for the SS mixture. [Reproduced with permission from 
ref. (Kim and Yamamoto, 2000).] 



0.5 Kinetic Lengths at Amorphous Interfaces 

If your goal is to clarify a phenomenon in the homogeneous liquid, the inclusion of an 
interface is not usually a good idea. The problem is that the interface will typically 
perturb the liquid structure significantly and, thus, so alter the phenomenon from that 
found in the bulk as to obscure any connection with the homogeneous situation. There 
is, of course, considerable interest in exactly such perturbed situations in the context 
of glass transitions actually taking place in confined geometries. This subject has been 
reviewed by a number of authors (Binder and Kob, 2005; Baschnagel and Varnick, 
2005). Within the artificial world of simulations, however, it is possible to contrive a 
surface that is structurally neutral by simply freezing the motions of some portion of 
the liquid. Such walls are an example of the imposition of a kinetic constraint. The 
interesting question is then to determine the length scale over which the influence of 
such a constraint is propagated into the unconstrained liquid. The idea of using an 
interface to establish the kinetic length scale was described (Butler and Harrowell, 
1991) in 1991 in the context of a lattice model of glassy kinetics. 

Scheidler et al (Schcidler et al, 2000; Scheidlcr et al, 2004) have reported simu- 
lation studies of the BMLJ1 mixture adjacent to a rough wall made up of the frozen 
liquid. These workers have considered both a cylindrical pore (Schcidler et at, 2000) 
and a planar wall (Schcidler et al., 2004). While the idea of a frozen liquid wall is 
simply sketched, its implementation takes some care. In ref. (Scheidler et al., 2004) 
the temperature of the liquid used to produce the frozen walls was adjusted to mini- 
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mize any structural perturbation and an additional repulsive potential was included to 
prevent particle penetration into the wall. Schcidlcr et al (Schcidlcr et al., 2004) fitted 
the relaxation time T q (z) (associated with the decay of a self intermediate relaxation 
time F s (q, z, t) at a distance z from the wall) for the BMLJ1 mixture to three different 
functions of the normal distance z from the planar wall. The best fit was obtained 
with the following expression, 

T q (z) 



In 



A(T) exp 



(0.16) 



UT)_ 

The kinetic length £ D increased by a factor of 3 as the temperature dropped to 
T = 0.5. The increase in £ Q with temperature was fitted with a simple Arrhenius form. 
The temperature dependence of £ appears to be quite different from the (T — T c ) 1 
dependence reported for the length of the dynamic heterogeneities in the bulk for the 
same model (Donati et al., 19996). Obtaining values of the surface kinetic length at 
temperatures closer to T c have been frustrated by poor statistics. 

Where Scheidler et al (Scheidler et at, 2000; Schcidlcr et al., 2004) have described 
how effective a frozen liquid is in imparting its immobility to an adjacent mobile liquid, 
Cavagna et al (Cavagna et al., 2007; Biroli et al., 2008) have examined how the frozen 
liquid can actually constrain the configurations available to the mobile liquid. The 
use of amorphous boundaries to establish an equilibrium length scale associated with 
structural correlations was described by Bouchaud and Biroli (Bouchaud and Biroli, 
2004) and Montanari and Semerjian (Montanari and Semerjian, 2006). Instead of a 
frozen half plane, Cavagna et al (Cavagna et al., 2007) have considered a spherical shell 
of frozen SS liquid, enclosing the mobile liquid in a volume of radius R. We emphasise 
that this study does not deal with the kinetics of the confined liquid but, instead, uses 
an amorphous wall to extract a static length scale. It is included in this review of 
kinetic lengths, in part, because of the obvious methodological parallels with the work 
of Scheidler at al (Scheidler et al., 2004). There is also a general expectation that kinetic 
length scales derive from underlying static length scales. Mark Ediger (Ediger, 2000) 
expressed the situation with admirable delicacy, "At present, it is an article of faith 
that something in the structure is responsible for dynamics that can vary by orders of 
magnitude from one region of the sample to another at T s " . Recent work (Cammarota 
et al., 101) has suggested that the origin of the growing kinetic length scale might lie 
in a separation of phases characterised by distinct degrees of the amorphous order. 
In ref. (Biroli et al, 2008), the authors determine, using an accelerated Monte Carlo 
algorithm, the degree of overlap q c of the ensemble of configuration at the centre of 
the sphere as a function of the radius R. They found that the overlap decayed to the 
random value q as 

q c (R)-q = ncxp[-(r/tf} (0.17) 

with both the length £ and the exponent £ increasing with supercooling - at T = 0.482 
the length £ = 0.617 and the exponent £ ~ 1 and both increase on cooling so that at 
T = 0.203, £ = 3.82 and C = 4.00. 

An earlier study of a liquid confined within a frozen liquid shell was presented by 
Sim et al (Sim et al., 1998; Sim et al., 1999). These authors studied a single component 
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Lennard- Jones liquid in 2D where the crystallization was frustrated by the disorder 
of the frozen wall. They did not present any systematic results associated with a 
kinetic length. However, the model is of potential interest due to the simplicity of the 
confined liquid and the possibility of an explicit counting of allowed configurations. 
The model represents an extension of the classic problem of packing disks in finite 
containers (Desmond and Weeks, 2009). Sim et al (Sim et al., 1999) reported some 
specific examples of collective reorganization events. 

The geometry of the pinned particles - planes, tubes or spherical cavities - reflect the 
particular problem that inspired the authors. In the case of Scheidler et al (Scheidler 
et al., 2004), this was the influence of confinement on glassy dynamics, while for 
Cavagna et al (Cavagna et al., 2007) it was to test the idea of droplet excitations as 
conceived in the mosaic theory (Kirkpatrick and Wolyncs, 1987; Kirkpatrick et al., 
1989; Xia and Wolynes, 2001; Lubchenko and Wolynes, 2003; Bouchaud and Biroli, 
2004). What about a spatial distribution of pinned particles that does not impose 
any particular correlation? Supercooled liquids subject to random pinnings have been 
studied by a number of groups (Kim, 2003; Lin et al., 2006). In 2003 Kim (Kim, 2003) 
reported the effect of pinning the positions of a fixed number Nd of randomly selected 
particles in the SS mixture. He found that the relaxation time r a scaled with Nd and 
T as 

r Q (T,iV d )(xcxp[iV d /T-] (0.18) 

with v = 3.7. In trying to extract a kinetic length scale from this calculation, Kim 
resorted to the following argument. There is one length scale imposed by the defects, 
i.e. 

&(T) cx (V/N d ) 1/3 (0.19) 

and there is another, the intrinsic kinetic length £(T), the one we are actually interested 
in. If one assumes that a) the intrinsic kinetic length is unperturbed by the pinned 
particles and b) that at the glass transition (defined as r a equalling some big number) 
\ d {T) = £(T), then it follows that 

£(T) oc T-" /3 (0.20) 

The argument is awkward. In particular, assumption b) above neglects the crossover 
to simple unpinned behaviour when £d(T) > £,(T). The methodology, however, has 
potential as a general tool for exploring kinetic and static correlations through the 
imposition of dilute random pinnings. 

0.6 Kinetic Lengths from Cross-Over Behaviour 

We shall complete our survey of length scales with those that arise most directly 
from the dynamical processes of interest and, therefore, perhaps represent the most 
pressing demand for our attention. The length scales described in this Section take 
as their starting point the existence of some sort of length dependent crossover in a 
physical property explicitly associated with particle motion. The existence of dynamic 
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heterogeneities leads, not surprisingly, to a range of physical behaviour that deviates 
from that expected of a uniform system. As one probes the behaviour over length 
scales larger than that of the heterogeneity, the 'classical' behaviour is recovered and, 
accordingly, a length scale can be associated with this crossover to homogeneity. 

Self (tracer) diffusion is one such phenomenon that exhibits a crossover that is 
manifest in changes in the time dependent displacement probability distribution G(r, t) 
(the van Hove distribution function). The story around G(r,t) in glass forming liq- 
uids is quite rich and so we shall take a moment to sketch out some main points, 
at least to clarify some terminology. In 1990, Odagaki and Hiwatari (Odagaki and 
Hiwatari, 1990) noted that, at fixed times, G(r, t) underwent a transition from Gaus- 
sian to non-Gaussian as a liquid was supercooled. Hurley and Harrowell (Hurley and 
Harrowcll, 1996) pointed out that this was an expected consequence of the increase in 
dynamic heterogeneity; specifically, the presence of persistent kinetic subpopulations. 
Dynamic heterogeneities have also been invoked to resolve another puzzle involving 
self diffusion. Sillescu and coworkers (Fujara et al., 1992) had shown experimentally 
that, on supercooling, fragile liquids exhibited a breakdown in the scaling between the 
diffusion constant D and both the shear viscosity (the Stokcs-Einstein relation) and 
the rotational diffusion constant (the Debye expression). This phenomenon is actu- 
ally evident just among the different length scales of the self intermediate scattering 
function F s (q, t), the Fourier transform of G(r, t). The temperature dependence of the 
relaxation time of F s (q,t) exhibits, on supercooling, an increasing anomalous depen- 
dence on q, with the small q time scale behaving like D (as it must) and the large 
q relaxation exhibiting a temperature dependence similar to that of the shear viscos- 
ity (Pcrcra and Harrowcll, 1998; Chandler et al., 2006; Chaudhuri et al, 2007). This 
general loss of a single time scale on supercooling is referred to as 'decoupling'. A 
number of groups (Cicerone and Ediger, 1996; Bcrthicr, 2004; Bcrthier et al., 20056) 
have reached the following consensus regarding the origin of this decoupling. The idea 
is that different transport properties correspond to different moments of the distribu- 
tion of microscopic times, so their decoupling at a particular wave vector is associated 
with the growth of dynamical heterogeneity (as manifest in the broadening of the 
distribution of microscopic times) over the corresponding length scale. Chaudhuri et 
al (Chaudhuri et al, 2007) have pointed out that the presence of an exponential tail in 
the van Hove function is a signature of the presence of slow and fast particles and can 
account for the decoupling of diffusion and structural relaxation. The continuous time 
random walk model they propose has been further quantified by Hedges et al (Hedges 
et al., 2007) who demonstrated that the ratio of the persistence time over the exchange 
time (the times for the first move and between subsequent moves, respectively) grows 
rapidly in the supercooled liquid. 

The explanation of decoupling provided in refs. (Chaudhuri et al, 2007; Hedges 
et al., 2007) requires that the different transport processes can be related to a common 
distribution of microscopic times. While this condition is met for the relaxation of 
F s (q, t) at different q, it is not clear that it is met by the Stokes-Einstein breakdown 
itself since diffusion and viscosity correspond to quite different physical processes and, 
therefore, are associated with quite distinct distributions of microscopic times. Other 
approaches (Hodgdon and Stillingcr, 1993) to the decoupling of diffusion and viscosity 
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Fig. 0.6 The time dependence of the mean square displacement for the A particles in the 
BMLJ1 model for T = 1.0, 0.8, 0.6, 0.55, 0.50, 0.47 and 0.45 listed from left to right. The 
symbols are placed at different characteristic times. Squares: the time at which the standard 
non-Gaussian parameter reaches the maximum value. Triangles: the a relaxation time r a . 
Circles: the onset time for Fickian diffusion. [Reproduced with permission from ref. (Szamel 
and Flenner, 2006).] 

avoid this particular criticism by retaining explicit coupling between local mobility 
fluctuations and stress relaxation. 

Over a long enough time, particles will sample a sufficient number of kinetic envi- 
ronments so as to recover standard or Fickian diffusion. Szamel and Flenner (Szamel 
and Flenner, 2006) have determined the time tf over which a particle must, on av- 
erage, move before G(r, t) becomes Gaussian for the BMLJ1 mixture simulated using 
Brownian dynamics. Their results are shown in Fig. 0.6. They observe that tf is larger 
than the relaxation time T a and that this difference increases as the temperature is 
lowered. A length scale can be obtained from Fig. 0.6 simply by reading off the mean 
square displacement at t = tf- This length I* ranges up to 2.5 small particle diam- 
eters and corresponds to the distance that a particle must on average move before 
exhibiting Fickian diffusion. The value of I* from ref. (Szamel and Flenner, 2006) is 
roughly half that predicted by the expression for the crossover length due to Bcrthicr 
et al (Bcrthicr et at, 2005c), I* cx \fDr a , but, at large supercoolings, the two lengths 
exhibit a similar temperature dependence. 

Stariolo and Fabricius (Stariolo and Fabricius, 2006), in a study of the BMLJ1 
mixture, reported the appearance of a new crossover in the self intermediate scattering 
function at around r Q , in addition to the crossover to Fickian diffusion at longer times. 
The authors associated the earlier crossover with the length scale of the dynamic 
heterogeneities as measured by %4. A link between the length scale of the crossover to 
Fickian behaviour and the length scale of the dynamic heterogeneities as obtained from 
the 4 point correlations has been explored by Bcrthicr (Bcrthicr, 2004). Studying the 
BMLJ1 model, he demonstrated that a suitably normalised product of the q-dependent 
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relaxation time and the diffusion constant from a wide range of temperatures could be 
collapsed onto a single master curve when the wavevector q was scaled by the length 
scale of dynamic heterogeneities. 

Is the crossover to Fickian behaviour governed by a time scale tf or a length scale 
1*1 Given the transient character of dynamic heterogeneities, the answer appears to 
be a matter of taste. A crossover involving an unequivocally static distribution of het- 
erogeneities has been studied by Barrat and coworkers (Leonforte et ai, 2005). These 
authors demonstrated that the continuum elastic description of a disordered polydis- 
perse mixture of Lennard- Jones particles at zero temperature broke down for length 
scales less than ~ 23 particle diameters. Over length scales smaller than this crossover 
length, the material exhibited non-affine response to an applied strain and this length 
gives the lower wavelength bound for the applicability of classical eigenvectors. 

The connection between non-affine displacements and the approach to the glass 
transition has been made explicit in a study by Mosayebi et al (Mosaycbi et ai, 2010). 
The local potential energy minima of the BMLJ1 model have been collected from MD 
trajectories as a range of temperatures. These inherent structures are subjected to a 
static shear deformation at T = and the spatial distribution of the resulting non- 
affine displacements calculated. These authors find that a characteristic length scale 
of the non-affine field grows as the temperature from which the inherent structures 
are obtained is lowered. 

While the complexities of normal modes of disordered materials has not tradi- 
tionally been associated with dynamic heterogeneities, a growing body of evidence 
suggests that the two phenomena are correlated (Schober et ai, 1993; Oligschleger 
and Schober, 1999; Brito and Wyart, 2006; Brito and Wyart, 2007; Widmer-Cooper 
et al., 2008; Widmer-Cooper et ai, 2009). It is possible that the large, but finite, length 
scales identified in ref. (Leonforte et at, 2005) represent a useful T = limit for the 
length scale of dynamic heterogeneities. 

So far in this review we have made no mention of a length scale associated specifi- 
cally with stress relaxation. Certainly, the overwhelming emphasis of simulation studies 
of the glass transition has been on single particle dynamics, in spite of the central role 
of viscosity in defining the glass transition. There are, however, a number of interesting 
papers on growing length scales associated with transverse momentum fluctuations, the 
kernel of the Green-Kubo expression for shear viscosity. In 1995, Mountain (Mountain, 
1995) used the transverse momentum autocorrelation function to obtain the disper- 
sion curve for the supercooled SS mixture. He identified a growing length scale as the 
longest wavelength associated with propagating transverse modes. This wavelength, 
obtained by extrapolating the dispersion curve for the supercooled liquid to the point 
where the mode frequency vanished, showed a strong increase at large supercoolings. 
At the temperature at which Yamamoto and Onuki (Yamamoto and Onuki, 1998a) 
found a kinetic length of ~ 10, Mountain reports a length of ~ 55 ( a value roughly 
three times the length of his own simulation box). Ahluwalia and Das (Ahluwalia and 
Das, 1998), using a mode coupling theory, have argued that the length identified by 
Mountain will diverge as T g is approached from above. There is a serious difficulty in 
associating the length described by Mountain and Das with the kinetic length scales 
that are the subject of this review. As pointed out by Hiwatari and Miyagawa (Hiwatari 
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Fig. 0.7 The wavevector-dependent viscosity for the BMLJ1 mixture divided by the k = 
value, vs the reduced wavevector for 200 K, 160 K, 120 K, 100 K, 80 K and 60 K. The value 
of f, the kinetic length, is obtained from the fitting function in Eq. 0.21. [Reproduced with 
permission from ref. (Kim and Keyes, 2005).] 

and Miyagawa, 1990), a straightforward application of viscoelastic theory results in the 
prediction that the longest wavelength associated with transverse propagation is pro- 
portional to rj, the shear viscosity. This relation simply reflects the condition that, for 
a mode to propagate, its frequency must exceed 1 / 77 , the value set by the dissipation. 
The growing length described by Mountain (Mountain, 1995) is a direct consequence 
of the growing relaxation time and is quite independent of any microscopic correlation 
length associated with the physical origin of this growing relaxation time. 

An alternate and more informative treatment of the transverse momentum corre- 
lation function has been presented by Kim and Keyes (Kim and Keyes, 2005). These 
workers have calculated the time integral of the k-dependent transverse momentum 
autocorrelation function. This quantity, essentially the zero frequency component of 
the correlation function, is the wavevector-dependent shear viscosity r)(k). The au- 
thors (Kim and Keyes, 2005) found that r)(k), calculated for the BMLJ1 model, could 
be fitted with the following functional form, 



1 + atanh{k£(T)) 



(0.21) 



Consistent with Eq. 0.21, r](k)/r](0) plotted against k£(T) collapses the data from 
different temperatures onto a single master curve (as shown in Fig. 0.7) that decays 
with increasing wavevector. The length scale £ was found increase from 0.13 to 1.62 
(taa on cooling. Kim and Keyes (Kim and Keyes, 2005) go on to argue that the 
increase in £ can be directly linked to the breakdown in the Stokes-Einstein scaling 
through the use of a mode coupling expression due to Keyes and Oppcnhcim (Keyes 
and Oppcnhiem, 1973). 
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Furukawa and Tanaka (Furukawa and Tanaka, 2009), studying the SS model, have 
extended the Kim-Keyes observations in two ways. First, they have shown that the 
increasing wavevector-dcpcndcncc of r)(k) is due entirely to the transverse component 
of the momentum flux and, second, they have established (again, for the SS mixture) 
that the length extracted from rj(k) has the same magnitude and temperature depen- 
dence as X4 obtained, as described in Section 0.3, from fluctuations in the structural 
relaxation function. Puscasu et al (Puscasu et al., 2010a; Puscasu et al., 20106) have 
simulated the wavevector dependence of the velocity kernel for a range of liquids in- 
cluding diatomic molecules and short chain alkanes and, for the latter model, report 
the growth of a very large length scale as the supercooling is increased. 

The decay of r)(k) with wavevector k reflects the decrease in dissipation as the 
wavelength shortens. As harmonic solids show no dissipation, it is tempting to as- 
sociate the growing length identified in refs. (Kim and Keyes, 2005; Furukawa and 
Tanaka, 2009) as being associated with the characteristic length scale of such solid- 
like domains. The (rough) superposition of the normalised r/(fc)'s from both high and 
low temperatures, as shown in Fig. 0.7, suggests that, in the supercooled liquid, stress 
relaxes similarly to that in the high temperature liquids except that the elementary 
objects are now rigid clusters with a linear dimension £ instead of the atomic compo- 
nents. In 1989, Ladd and Alder (Ladd and Alder, 1989) described the stretched tail of 
the shear stress autocorrelation function in hard sphere liquids near freezing. (Their 
evocative label - the 'molasses' tail - does not seem to have caught on.) Subsequently, 
Isobe and Alder (Isobe and Alder, 2009; Isobe and Alder, 2010) have argued that the 
long time relaxation of the shear stress is dominated by the life time of rigid clusters 
in the liquid. 

0.7 What Lengths Influence Relaxation? 

Does a kinetic length provide the unified description of cooperative dynamics for which 
it was intended? We have seen that the relaxation time of a glass forming system can 
be increased by either decreasing the temperature or decreasing the system size below 
some threshold value. Karmakar et al (Karmakar et al., 2009) have examined the 
dependence of both the relaxation time r and the susceptibility Xi{t*) {Xi m their 
notation) of the BMLJ1 mixture as a function of temperature and number of particles. 
Their results are shown in Figs. 0.8 and 0.9. Fixing the system size at a large value, 
say N = 1000, we see that Xa an( l r both increase as T decreases, similar to behaviour 
already described in Section 0.3. If, however, we hold T fixed, the variation of Xa an< ^ 
r with respect to N have opposite signs. The authors note that this result is contrary 
to the expectations of finite scale scaling. It indicates that xf does not contain all of 
the information about the collective processes in the liquid necessary to establish the 
relaxation time. 

The puzzling observations of Karmakar et al (Karmakar et al., 2009) had been fore- 
shadowed by earlier work. Kim and Yamamoto (Kim and Yamamoto, 2000) demon- 
strated that the increase in the relaxation time they observed for the small system 
was accompanied, not unexpectedly, by a truncation of the size of mobile clusters. 
The point here is that the more extended mobile regions may represent greater mo- 
bility for the system, not less. Jackie and coworkers (Frobose et al., 2000; Jackie and 
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Fig. 0.8 Peak height of the dynamic susceptibility, Xi (T,n) for the BMLJ1 model plotted 
as a function of system size N for different temperatures. For each temperature, Xi{T,ri) 
increases with system size, and saturates for large system sizes. Xi{T,n) also increases as 
the temperature is lowered. Insert: Xa{T, n) plotted as a function of time. In the main plot, 
temperature increases from the top curve to the bottom. In the insert, temperature decreases 
moving from the left curve to the right. [Reproduced with permission from ref. (Karmakar 
et al, 2009).] 

Kawai, 2001) demonstrated that, even in systems exhibiting dynamic heterogeneities, 
the diffusion of particles still involved significant coupling to extended visco-elastic 
flows (i.e. collective motions in which particles retain their neighbours) dominated by 
the longest wavelength transverse modes in the system. Decreasing the system size 
removes the longer wavelength modes and, possibly, stiffens the surrounding medium, 
leading to a decrease in mobility. Such a scenario would be expected to exhibit a 
very different system size dependence to the way size influenced the dynamic hetero- 
geneities themselves. The spatial correlations in small displacements, as opposed to 
the large ones that typically define 'mobile particles, are an important component of 
the length that characterises relaxation, both incoherent (e.g. self diffusion) and coher- 
ent (e.g. stress relaxation). In a recent study of structural relaxation (Widmcr-Coopcr 
et al., 2009), it was shown that, of the particle movements that have contributed irre- 
versibly to relaxation, 60% (at the lowest temperature studied) could be categorized 
as strain-like, meaning that they involved the loss of no more than one of the initial 
neighbours. 

0.8 Conclusions 

The measures of a kinetic length scale reviewed in this chapter have succeeded in a 
number of things. They have confirmed, by the act of measurement, the existence of 
spatial heterogeneity of the kinetics and the coarsening of this distribution on cooling. 
This descriptive success is an important milestone. It allows for the comparison be- 
tween different glass formers, between the spatial character of the dynamics and that 
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Fig. 0.9 Relaxation times as a function of temperature and system size in the BMLJ1 
mixture. For the smallest temperature, r(T, N) increases by approximately a decade from 
the largest to the smallest system size. Temperature increases from the top curve to the 
bottom. [Reproduced with permission from ref. (Karmakar et ai, 2009).] 

of the static properties of the supercooled liquid and between various theoretical treat- 
ments of the glass transition. The repeated observation of power law relations between 
a kinetic length and a relaxation time is significant. Such behaviour is suggestive of 
relaxation processes governed by the fluctuations in domain size. An alternative, in 
which the length corresponded to the dimension of an object involved in an activated 
process, would be expected to exhibit a different relation (Cavagna, 2009), i.e. 



The power law exponents have been found to lie, roughly, between 2 and 4. If we 
think of the glass transition as being the temperature at which a relaxation time as 
increased by a factor of, let us say, 10 10 , then these power laws would require the kinetic 
length at this glass transition to have increased by a factor of between 10 2 and 10 5 . 
The reports of temperature dependences for the kinetic length with singularities at the 
mode coupling temperature are generally regarded as a signature, not of a divergence, 
but of a cross-over to some alternative relaxation mechanism. To date there is no 
evidence for the large increases in kinetic lengths suggested by the observed power 
law. The conclusion is then that the time-length relationship extracted over a limited 
range of temperatures in simulations does not continue to hold at lower temperatures. 
The slow down in the growth of the kinetic length on cooling has been noted in 
experiments (Dalle-Ferrier et ai, 2007; Brambilla et ai, 2009). In simulations of the 
BMLJ1 mixture, Berthier et al (Berthier et ai, 20076) reported that for T < 0.47, the 
growth of the dynamic susceptibility with respect to r Q becomes much slower than that 
observed at higher temperatures, "perhaps logarithmically slow" . This last comment 
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from (Berthier et al., 2007 6) is a reference to the prediction of a logarithmic relation 
£ ~ (lnr Q )^ , such as expressed in Eq. 0.22, from theories that invoke an activated 
process (Lubchcnko and Wolynes, 2003; Bouchaud and Biroli, 2004). 

The kinetic lengths reviewed in Sections 0.2 and 0.3 all provide useful and, es- 
sentially similar, measures of the dynamic heterogeneities. An important question is, 
however, how well can they account for the length scales implicit in finite size effects, 
interfacial correlations and crossover behaviour? Furukawa and Tanaka (Furukawa 
and Tanaka, 2009) have reported that the characteristic length obtained from the 
q-dependence of the transverse momentum fluctuations exhibits a similar size and 
temperature dependence to Xi{t*)- Berthier's (Berthier, 2004) demonstration that the 
Fickian crossover can be scaled by the length scale from S±{q,t) establishes a similar 
link. On the other hand, the increase in time scale due to the reduction in system 
size points to the role of correlations not included in the dynamic susceptibility. There 
remains an open problem to establish a clear identification of an explicit kinetic length 
(i.e. analogous to the lengths defined in Sections 0.2 and 0.3) with an implicit length 
scale such as demonstrated in Sections 0.4, 0.5 and 0.6. 

As already touched upon, there is a gap between description and mechanism, one 
that only becomes truly evident now that the problem of description has been largely 
solved. To date, it is not clear that dynamic heterogeneities have clarified the mecha- 
nisms of relaxation. Part of the challenge in addressing the issue of mechanism certainly 
lies in refining what it is we actually want explained. The growth of the activation en- 
ergy on cooling fragile liquids is generally associated with the growth of the number 
of elementary processes that must occur in series to achieve relaxation. There is cer- 
tainly clear evidence, through the finite size results, amorphous interface studies and 
the observation of the characteristic length in the transverse momentum fluctuations, 
that the mechanisms responsible for relaxation of particle positions and stress do have 
characteristic lengths and that these grow as the temperature decreases. Missing is 
a description of those elementary process by which the observed length scales are 
generated. 

There remain a number of interesting and open challenges that arise directly from 
what we have already learnt from MD simulations. We shall conclude by listing a 
few. First, there is evidence that the number of mobile particles (however one might 
define them) in a supercooled liquid undergoes substantial fluctuations in time. The 
democratic particle approach, introduced by Appignanesi and coworkers (Appignanesi 
et al., 2006a; Appignanesi et al., 20066), provides one explicit measure of the inter- 
mittent fluctuations in the number of particles involved in large displacements. From 
experiments on granular systems (Candalier et al., 2009) and, subsequently, MD simu- 
lations (Candalier et al., 2010), Dauchot and coworkers have identified the intermittent 
appearance of spatially correlated bursts of enhanced mobility among particles. Chris- 
tened 'avalanches', this intermittent behaviour appears to increase as the glass transi- 
tion is approached. Clarification of the significance of these fluctuations is important. 
Is there a new length scale (or a hierarchy of new lengths) describing the volume over 
which the avalanche occurs? Is relaxation at low temperatures increasingly dominated 
by intermittent bursts of activity and, if so, what is happening during the quiescent 
intervals that leads to initiating a mobility event? 
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A second challenge is to make direct connection between dynamic heterogeneities 
and coherent processes like shear stress relaxation. The growing length scale associated 
with the transverse momentum fluctuations (Kim and Keyes, 2005; Furukawa and 
Tanaka, 2009; Puscasu et ai, 2010&) provides a starting point. Along with particle 
mobility, dissipation is also becoming increasingly heterogeneously distributed as the 
temperature drops. To understand such phenomena we need to understand how to 
construct a description of visco-elastic behaviour where the dissipation can be localised 
even as the elastic behaviour becomes more extensive with the approach to the glass 
state. 

Finally, it remains a fundamental tenant that the ultimate origin of a kinetic length 
scale lies in length scales associated with structural correlations. Understanding this 
link between structure and dynamics is a problem that, currently, only simulations can 
address. There are no shortage of aspirants for the missing link: the spatial distribution 
of localised soft modes (Schober et ai, 1993; Oligschlcgcr and Schobcr, 1999; Brito and 
Wyart, 2006; Brito and Wyart, 2007; Widmer-Cooper et ai, 2008; Widmer-Coopcr 
et at, 2009), the mosaic length scale (Cavagna et al, 2007; Biroli et ai, 2008), the 
physical extent of clusters of locally preferred structures (Coslovich and Pastorc, 2009; 
Tanaka et ai, 2010; Lerner et ai, 2009; Mossa and Tarjus, 2006; Pedersen et ai, 
2010) to name some. Should the relevant structural length scales ever be successfully 
unearthed it is quite possible that the description of dynamic heterogeneities may be 
rendered irrelevant, replaced by the more convenient and illuminating account provided 
by the relevant structure. That dynamic heterogeneitiy may prove the agent of its own 
demise (as a descriptor of coopcrativity, that is) is a real possibility. Until then, the 
spatial distribution of dynamics remains our most general and concrete description of 
the complex dynamics associated with a liquid's passage to rigidity in the glass state. 
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